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ABSTRACT 

The  theory  of  small  superposed  deformations  for 
isotropic  incompressible  elastic  materials  is  used  (i)  to  obtain 
necessary  restrictions  on  the  form  of  the  strain -energy  function 
by  requiring  that  the  speed  of  propagation  be  real  for  waves  that 
pass  through  a finitely  deformed  body  of  material  (i.e.,  Hadamard 
stability  criterion) , and  (ii)  to  determine  critical  loading  con- 
ditions for  a thick  rectangular  plate  under  which  bifurcation 


solutions  (i.e.,  adjacent  equilibrium  positions)  can  exist.  The 
possibility  of  bifurcation  under  tensile  loading,  when  one  pair  of 
faces  of  a plate  are  force  free,  is  precluded  by  further  material 
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1 . Introduction 

The  main  purpose  of  this  paper  is  to  bring  together  results 
obtained  during  the  past  few  years  on  two  basic  problems  in  finite 
elasticity  theory.  One  of  these  is  the  question  of  what  restric- 
tions  must  be  imposed  on  a constitutive  assumption  in  order  to  en- 
sure physically  realistic  behavior  of  the  material  which  it  claims 
to  describe.  The  other  problem  deals  with  the  description  of  the 
conditions  under  which  adjacent  equilibrium  states  can  appear  as 
small  perturbations  on  an  equilibrium  state  of  finite  deformation 
of  a thick  plate.  Most  of  the  results  appear  in  references  [ 1 ] - [ 5 ] . 

The  setting  for  both  problems  is  the  theory  of  small  deforma- 
tions superposed  on  states  of  pure  homogeneous  deformation  of  bodies 
made  of  an  isotropic,  incompressible  elastic  material  [6].  The 
relevant  equations  are  beiefly  derived  in  §2. 

The  approach  taken  for  the  first  problem  is  to  investigate  the 
implications  that  stem  from  the  Hadamard  stability  criterion,  which 
requires  that  the  speed  of  propagation  be  real  for  plane  waves  that 
pass  through  a deformed  body  of  material.  It  was  shown  by  Hayes 
and  Rivlin  [7]  that  this  criterion  must  be  imposed  in  order  to  en- 
sure that  the  body  can  actually  exist  in  the  finitely  deformed  state 
considered.  Accordingly,  we  adopt  the  resulting  implications  as 
minimal  material  stability  conditions . 

The  equations  that  describe  wave  propagation  are  derived  in 
§3,  and  further  discussed  in  §4,  for  a general  material.  The  main 
result  is  that  certain  inequalities  [1]  {see  eqn.(4.1)  below}  in- 
volving derivatives  of  the  strain -energy  function  are  necessary  and 
sufficient  for  the  Hadamard  criterion  to  be  satisfied  for  waves 
whose  directions  of  propagation  lie  in  a principal  plane  of  strain. 
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An  open  question  remains  as  to  whether  these  are  also  necessary 
and  sufficient  when  a completely  general  propagation  direction  is 
considered.  It  is  seen,  however,  that  the  answer  to  this  is  affir- 
mative when  two  principal  extension  ratios  are  equal,  or  in  §5, 
when  the  constitutive  assumption  is  specialized  [3], 

The  determination  of  critical  states  of  finite  deformation, 
at  which  non-trivial  adjacent  equilibrium  states  first  become  pos- 
sible, has  been  the  subject  of  numerous  investigations  over  the 
past  twenty  years.  A partial  list  of  these  is  contained  in  [4]. 

The  geometry  considered  is  a rectangular  plate,  one  pair  of 
whose  faces  are  force  free  and  whose  other  two  pairs  are  loaded  by 
uniformly  distributed  normal  forces  when  the  body  is  held  in  its 
finitely  deformed  state.  The  forces  acting  on  one  opposite  pair  of 
faces  are  regarded  as  passive,  merely  holding  the  length  of  the 
plate  fixed  in  that  direction.  The  forces  acting  on  the  other  pair 
are  regarded  as  being  the  "cause"  of  the  bifurcation,  which  is 
assumed  to  appear  as  a small  plane  static  deformation  parallel  to 
the  faces  bearing  the  "passive"  forces. 

The  equations  governing  such  bifurcations  are  derived  in  §6, 
and  it  is  found  that  further  discussion  of  the  problem  can  be  con- 
veniently carried  out  if  two  cases  that  arise  naturally  are  separa- 
ted. This  separation  is  governed  by  the  particular  values  assumed 
by  the  one  relevant  material  parameter  that  appears  {see  (6.15) 
and  (6.18)}.  One  case  [4]  is  considered  in  §7  and  the  other  in  §8. 

When  the  causal  force  is  a thrust,  a single  asymptotic  formula 
[5]  {see  (8.17)},  relating  the  cross-sectional  geometry  of  a plate 
to  the  critical  condition  for  flexural  bifurcation,  is  found  to  be 
valid  for  all  materials.  This  result,  for  a "thin"  cross  section, 
covers  the  case  of  classical  Euler-type  buckling. 
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When  the  causal  force  is  tensile,  it  is  found  that  a bifur- 
cation can  occur  only  if  the  load  decreases  after  attaining  a maxi- 
mum. However,  this  behavior  rules  out  the  possibility  that  a plate 
of  material  can  be  held  in  stable  equilibrium  in  states  of  pure 
homogeneous  deformation  just  prior  to  that  required  for  the  bifur- 
cation, and  leads  to  the  imposition  of  further  restrictions  on  the 
strain-energy  function  beyond  those  discussed  in  §4. 

A major  question  that  remains  open  for  the  bifurcation  problem 
considered  here,  particularly  for  "thick"  cross  sections,  is  whether 
the  body  actually  goes  over  into  a mode  shape  that  is  assumed  at 
the  point  where  the  thrust  load  reaches  a critical  value.  Some 
experimental  evidence  [8]  throws  doubt  on  the  matter.  Attempts  to 
answer  this  analytically  by  employing  energetic  considerations  have, 
so  far,  been  unsuccessful. 

Recently  [9],  [10],  the  bifurcation  problem  has  been  exhaus- 
tively studied  for  a plate  of  incompressible  material  by  employing 
an  incrementally-linear , time-independent  constitutive  relationship. 
This  theory  includes  small  superposed  elastic  deformations  as  a 
special  case.  Accordingly,  an  alternative  description  of  the  res- 
ults in  § § 7 and  8 could  be  gleaned  from  those  in  [9]  and  [10]  by 
expressing  the  relevant  incremental  moduli  in  terms  of  the  strain- 
energy  function. 
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2 . Basic  Equations 

We  let  £,  X and  x denote  the  vector  positions,  relative 
to  a fixed  rectangular  Cartesian  frame  x,  of  a typical  particle 
of  the  body  in  an  undeformed  state,  in  a finitely  deformed  state, 
and  in  a state  at  time  t,  respectively,  and  write 

x = X + u (2.1) 

where  u is  the  displacement  associated  with  an  additional  small 
superposed  deformation.  We  consider  states  of  pure  homogeneous 
deformation  whose  principal  directions  are  parallel  to  the  coordi- 
nate axes  of  x and  which  are  completely  characterized  by  the 
(positive,  constant)  extension  ratios  ’ Then 

X = FC,  F « diag(A1,A2,A3)  . (2.2) 

The  components  of  E, , X,  x and  u,  referred  to  x,  are  denoted 
by  £ , and  u^,  respectively.  In  the  sequel,  all  ex- 

pressions involving  u (and  its  derivatives)  shall  be  systemati- 
cally linearized  in  u . 

The  Finger  strain  matrix  B is  given  by 

B = FFT  = F2+  fF  + FfT  , (2.3) 

where,  with  (2.1)  and  (2.2), 

!'?*!•  f - II  ui,J  • <2-4> 

(We  employ  the  notation  ,a  to  denote  9/3£a  * as  well  as  the 
usual  summation  convention  for  repeated  lower-case  subscripts.) 

The  invariants  I1  and  I2  of  B are  given  by 

. — ::  ...  '■  — 1 
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Xi  = tr  ? = V ^ ’ 


I2  = hi (tr  B)2  - tr  B2}  = Ig+  ig  , 


(2.5) 


where 


ll  = Xl+X2+X3  ’ 12  = X2X3  + X3X1  + X1X2  ’ 

(2.6) 

3 3 2 

= 2 / XAUA,A  » 12  = 2 / XA(Il‘XA)uA,A  * 

(2.7) 

A=1  ’ A=1 

The 

incompressibility  condition,  det  F = 1,  together  with 

(2.2) 

and 

(2.4),  yields 

x1x2x3  = 1 

(2.8) 

and 

3 

I/a  ua,a  = 0 • 

(2.9) 

A = 1 ’ 

The  strain-energy  density,  W,  depends  on  B through 

h 

and  i0  . Letting  II.  denote  the  components  of  the  Piola- 
Kirchhoff  stress  and  P a hydrostatic  pressure,  we  have 


where 


* ^P£ljkcc6YFj6Fky  = 3W/3Fia 

' H(VWFicT  ViBfjsV  ' 

is  the  alternating  symbol  and 


WA  = 9W/3Ia  , A=1 ,2  . 


(2.10) 


(2.11) 


Following  the  previous  notational  conventional,  we  write 


J 
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Scd  = nai  + %i  • P = P + P , (2.12) 

where  and  p are  the  stress  and  hysrostatic  pressure,  res- 

pectively, associated  with  a superposed  deformation.  Also,  we  ex- 
pand WA  in  a series  about  I1  and  ; thus, 

"a  ■ WA  * WAlll  * ■ A=1-2  (2-13> 

where  WA  and  WAg  denote  the  first  and  second  derivatives  of 
W with  respect  to  IA  , evaluated  at  IA  = IA  . The  substitution 
from  (2.12)  and  (2.13)  into  (2.10),  with  (2.2),  (2.4),  (2.6)  and 
(2.8),  yields 

(2.14) 


(2.15) 


(2.16) 


(No  summation  is  implied  for  repeated  upper-case  subscripts.) 


nAA  = 2XA{W1+(I1-VW2}  - P/XA  * 
nAB  = o , A + B , A,B  = 1,2,3 

and 

\a  =b^  {cab6ab  +(i5ab'1)dab+  eab}ub,b*  p/xa  » 

nAB  = CABUB, A + °ABUA,B  » A **  B » 
where  <5^  is  the  Kronecker  delta  and 

CAB  = 2{W1+(I1*XA"XB)W2}  » 

°AB  = ^P‘2XAXBW2^/XAXB  * 

EAB  = 2XAXB{W2+2[Wll+(2Il-XA-X^Wl2 
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We  neglect  all  body  forces.  Then,  bearing  in  mind  that  the 
A's  are  constants,  from  (2.1),  (2.2),  (2.12)  and  (2.14),  the 
equations  of  motion  are 


IT  . = PU-  , 

ai ,a  l ’ 


(2.17) 


together  with  the  requirement  that  P be  constant,  where  p is 
the  mass  density  and  a superposed  dot  denotes  differentiation  with 
respect  to  t . The  substitution  from  (2.15)  into  (2.17)  yields 


B = 1 


tCABUA,BB  + EABUB,BA) 


P’A/XA 


pu, 


A = 1,2,3 


(2.18) 


Let  v denote  the  (outward  drawn)  unit  vector  normal  to  a 
surface  of  the  body  in  t'  formed  state.  Then  the  traction 

T acting  on  this  surf  Lime  t , and  measured  per  unit 

undeformed  area,  has  components  given  by 


where 


T . = n . v = T . + t . , 
i ai  a l i 


T.  = n -v  , t.  = n .v 
l ai  a ’ l ai  a 


(2.19) 


(2.20) 


r 
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3.  The  Propagation  of  Plane  Waves  : The  Secular  Equation 

In  this  section  we  derive  the  secular  equation  that  governs 
the  propagation  of  small  amplitude  sinusoidal  plane  waves  through 
a block  of  material  that  is  otherwise  held  in  equilibrium  in  the 
finitely  deformed  state  characterized  by  . Letting  to 

denote  the  angular  frequency,  < the  wave  number  and  n the  (real) 
unit  vector  in  the  direction  of  propagation  of  such  a wave,  we 
employ  the  usual  complex  notation  and  write 

u = U exp{  i (icn*X-<jot) } , (3.1) 

. p 

where  U is  a constant  vector  and  i = -1  . 

The  substitution  from  (3.1)  into  (2.9),  with  (2.2),  yields 

n»U  = niUi  = 0 , (3.2) 

where  we  disregard  the  trivial  case  k = 0 . 

We  express  the  hydrostatic  pressure  p associated  with  the 
superposed  deformation  in  a form  that  is  consistent  with  (3.1); 
thus , 

p = i<q  exp{ t (<n*X-ut)  } , (3.3) 

where  q is  a constant.  With  the  notation 

a = k2/  (poo2)  , (3.4) 

the  substitution  from  (3.1)  and  (3.3)  into  (2.18)  yields 

^(Qijnj  - qni)  = Uj  , (3.5) 


where 


Qij  = QijkUk  ’ 


(3.6) 


10 


with 


QABC  ' CABXBnB5AC  + XAXCnCEAC6AB 


(3.7) 


To  eliminate  q from  (3.5)  we  multiply  each  term  of  the 
equation  by  esrinr  an<*  obtain,  with  (3.6), 

Qstut  = 0 * 

where 

^st  " £srinr  ^a(^ij  tnj  0it^ 


(3.8) 


(3.9) 


Now  since  nsQst  vanishes  identically,  we  note  that  the  three 
equations  in  (3.9)  are  not  linearly  independent.  Taking  s = 1 
and  2,  for  example,  we  obtain 


«ltUt  * 0 • «2tUt  ■ 0 > 


(3.10) 


which,  with  (5.2),  constitutes  a set  of  three  homogeneous  equa- 
tions for  U . A non-trivial  solution  for  U can  exist  only  if 


£ijkQliQ2jnk  = 0 ’ 


(3.11) 


The  substitution  from  (3.9)  into  (3.11),  together  with  some 
algebraic  manipulation,  yields  the  desired  secular  equation  in 
the  form  [3] 


acr  - Ba  + 1 = 0 , 


(3.12) 


where 


a = £ijknknsnt^nl^2si^3tj+n2^3si^ltj  + n3^lsi^2tj-*  ’ 

0 = ns^jsj-njnkQjsk)  ' 


(3.15) 


and  where  we  have  used  the  fact  that  n*n  = 1 . According  to 
(3.4),  the  quantity  pa  is  the  square  of  the  slowness  associated 
with  a wave  described  by  (3.1).  We  see  that  (3.12)  yields  two 
positive  roots  for  a if  and  only  if 

a > 0 , 3 > 0 , and  g2-4a  >.  0 . (3.14) 

Before  turning  to  a discussion  of  the  implications  of  (3.14), 
we  must  obtain  explicit  expressions  for  a and  g in  terms  of  W. 
To  this  end  we  introduce  the  notation 

KA  “ V XAW2  * MA  = 2(W11+2XAW12+  XAW22)  ’ A=1>2>3  (3.15) 

and  also 

= ^nBnc(^g'^^)  > A / B ^ C / A • (3.16) 

Then  employing  (2.16)  and  (3.7)  in  (3.13),  with  (2.8),  we  obtain  [3] 
a = 4 (X2n2+  ##)  {(X^K2K3+  ##)  + W]_(J1+  ##)} 

+ 4W2(\2x2Ji+  ##) 

+ 16n2n^n2(X^-X^)2(X3-X^)2(X^-X2)2(W11W22-w22)  , (2.17) 

3 = 2{K]_(X2n2  + X2n2)  + J±  + ##  } , 

where  the  symbol  ##  denotes  two  additional  terms  which  are 
obtained  from  the  ones  shown  by  a cyclic  permutation  of  the  sub- 
scripts 1,2,3  on  the  X's,  n's,  K's,  M's  and  J's.  From  (3.17), 
with  (2.8),  (3.15)  and  (3.16),  we  calculate  B2-4a  : 
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= 4W2{[AjnJ(A2-A2)2+  ##] 

+ 2[A2A2n2n2(A2-A2)(A2-A2)  + *#)} 

+ 8W2{J1[A|(A2-A2)(n^+n2) 

+ ^^I'^Hn^n2)]  + ##} 

+ 4{ [ J1+  #*]2 

- 16n2n2n2(A2-A2)2(A2-A2)2(A2-A2)2(W11W22-W22)}  . (3.18) 

2 

' interest  to  note,  in  general,  that  the  quantity  6 -4a 
etely  independent  of  the  value  assumed  by  . 
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4 . Stability  Restrictions  for  a General  Material 

The  set  of  inequalities  (3.14),  which  must  hold  for  all  n , 
is  the  Hadamard  stability  criterion.  In  turn,  with  (3.17)  and 
(3.18),  these  inequalities  lead  to  restrictions  on  the  functional 
form  of  W , which  we  shall  adopt  as  minimal  material  stability 
conditions . 

It  has  been  shown  [1],  [3]  that,  when  the  propagation  direc- 
tion n lies  in  a principal  plane  of  the  underlying  pure  homogene- 
ous deformation,  the  necessary  and  sufficient  conditions  for  (3.14) 
to  be  satisfied  are 

Ka  > 0 and  (I1-Xj-2X-1)MA  + Ka  > 0 , A = 1,2,3  , (4.1) 

where  the  K's  and  M’s  are  defined  by  (3.15).  A sketch  of  the 
proof  of  this  will  be  given  below. 

The  first  set  of  conditions  (4.1)  is  known  as  the  Baker- 
Ericksen  inequalities  [11].  The  second  set  was  apparently  first 
given  in  [1].  It  was  shown  in  [2]  that  (4.1)  can  also  be  derived 
by  considering  small  static  superposed  shearing  deformations  in 
principal  planes  and  requiring  the  corresponding  incremental  shear 
moduli  be  positive. 

It  is  not  known  whether  any  further  restrictions  on  W , 
beyond  those  expressed  by  (4.1),  stem  from  (3.14)  for  an  arbitrary 
direction  n . However,  it  is  relatively  easy  to  show  that  (4.1) 
is  sufficient  to  ensure  8 > 0 for  arbitrary  n . We  employ 
(2.6)2,  (2.8)  and  (3.16)  in  (3.17)2  and  write 
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e = 2K1{(A2n^-X3n^2  + n2  (X2n2  + X2n2) } + **  + 2b  , 

009  21  (4.2) 

b = (X2  + X3)^n2n-{(I1-Xj-2X-x)M1+  + **  . 

From  (4.1)  we  see  that  b >.  0 and,  whence,  3 > 0 . 

In  the  general  case  when  all  three  X's  are  distinct,  we  note 
that,  at  most,  only  two  of  the  three  Baker -Ericksen  inequalities 
are  independent.  For,  with  the  ordering  x^  > X^>  X^  , say,  from 
(3.15)^  we  have  that  K1  > 0 and  K3  > 0 implies  K.^  > 0 . If 
two  of  the  X's  are  equal,  say  X^  * X2  j1  X^  , then,  again,  at 
most  two  of  the  Baker-Ericksen  inequalities  are  independent  while 
the  second  set  in  (4.1)  reduces  to  K3  > 0 and  just  one  other  in- 
equality. Thus,  for  these  special  types  of  deformations,  (4.1) 
yields,  at  most,  three  independent  restrictions  on  the  form  of  W . 
Some  additional  comments  pertaining  to  this  case  are  made  at  the 
conclusion  of  this  section. 

For  propagation  in  a principal  direction,  say  n = (1,0,0)  , 
(3.17)  yields 

a = 4X^K2K3  , 3 » 2X2(K2+K,)  , (4.3) 

and  it  is  clear  that  S2-4a  > 0 . From  expressions  for  a and  3 
analogous  to  these  for  propagation  along  the  other  two  principal 
directions,  we  can  conclude  that  a > 0 and  3 > 0 if  and  only 
if  > 0 . Essentially  the  same  results  can  be  read  off  directly 

from  results  obtained  much  earlier  by  Ericksen  [12]. 

>Ve  turn  now  to  the  case  of  propagation  in  an  arbitrary  direc- 
tion in  a principal  plane.  For  definiteness,  suppose  that 
n = (n. ,n2,0)  . Then  from  (3.17)  we  obtain 

a = 4y1Y2  , S s 2(y1+Y2)  , 

^ — .... . m 


(4.4) 
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where 

Y1  = XiniK2  + X2n2Kl  ’ 

y2  = K3CX1n2-X2n^2  (4.5) 

+ n2n2(X1+X2)2{K3+(X1-X2)2M3}  . 

It  is  clear  from  (4.4)  that  ?2  • 4a  > 0 . With  the  validity  of 
(4 . 1) ^ already  established,  it  easily  follows  [1],  [2]  from  (4.4) 
and  (4.5)  that  a > 0 if  and  only  if 

C X x - X 2 ) 2M 3 + K3  > 0 , (4.6) 

which,  then,  is  also  sufficient  for  8 > 0 , as  noted  previously. 
The  form  of  the  conditions  (4.1)^  follows  from  (4.6),  with  (2.6)^, 
by  considering  analogous  results  for  propagation  in  the  other  two 
principal  planes. 

If  the  underlying  deformation  is  a uniform  two-dimensional 
extension  or  compression,  we  have  a situation  where  two  of  the 
X’s  are  equal.  For  example,  say 

X1  = Xg  = X , X 3 = 1/X2  . 

Then,  with  (3.15),  the  inequalities  (4.1)  become 

K > 0 , K3  > 0 , (X3-1)2M  + XUK  > 0 , 

where 

K = K1  = K2  = W1+  X2W2  , K3  = Wx  * X‘UW2  , 

M - Mx  - M2  = 2(WX1  + 2X2W12+  XUW22)  . 


(4.7) 


(4.8) 


(4.9) 
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The  derivatives  of  W that  appear  in  (4.9)  are  evaluated  for 

I = 2X2  + 1/X1*  , I = Xk  + 2/X2  , (4.10) 

which  is  a point  on  the  curve  bounding  the  domain  in  the  I^g-plane 
where  W is  defined.  (See  [13],  for  example). 

Now,  when  two  principal  extension  ratios  are  equal,  we  note 
that  any  direction  n is  parallel  to  a principal  plane,  and  it 
follows  that  the  satisfaction  of  (4.8)  is  necessary  and  sufficient 
that  (3.14)  be  satisfied  for  arbitrary  directions  of  propagation 
[3]. 
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5 . Results  for  Some  Special  Materials 

A stumbling  block  in  the  way  of  obtaining  the  full  implica- 
tions of  (3.14)  for  arbitrary  n is  the  presence  of  the  term 

p 

(•W11W22"W12^  t*ie  exPressi°n  for  <*  • Three  special  cases  where 

this  term  is  absent  have  been  studied  [3] . These  are  the  Mooney- 
Rivlin  material  (see,  also,  [12])  and  materials  which  depend  either 
on  I or  on  I2  . For  these  it  is  found  that  the  basic  restric- 
tions implied  by  (4.1)  are  necessary  and  sufficient  for  (3.14)  to 
be  met  for  arbitrary  n . Only  the  case  where  W depends  on  I 
shall  be  discussed  here. 

We  write  W = f (I  ) so  that,  from  (3.15)  , 

KA  - f’  , Ma  = 2f " , A - 1,2,3  (5.1) 

and,  with  (2.6)^,  the  conditions  (4.1)  become 

I 

f'  > 0 , 2(XA-XB)2f"  + f'  > 0 , AB  = 12,23,31  . (5.2) 

The  derivatives  of  f are  evaluated  at  I1  = Ix  . From  (3.17), 

(3.18),  with  (5.1),  we  obtain 

a = 4(X2n2+  ##)f’{(X2n2+  ##)f' 

+ 2f"[(X2-X2)2n2n2+  ##]}  , j 

8 = 4{(X2n2+  ##) f'  + f" [ (X2-X2) 2n2n2+  ##]}  , (5.3) 


32-4a  = 1 6 { (X2-X2)  2n2n2+  ##}2(f")2  . 

It  is  clear  that  B2-4a  >.  0 , for  arbitrary  n , regardless  of 
the  form  of  f . Recalling  that  > 0 is  necessary  and  sufficient 
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that  the  speed  of  waves  propagating  in  a principle  direction  be 
real  {cf.(4.3)}  , it  follows  from  (5.1)1  that  (5.2)^  must  be  satis- 
fied. 

Now  if  f"  > 0 , we  see  that  (5.2)2  is  satisfied  and  also 

a > 0 . Therefore,  we  need  consider  only  such  f where  f"  < 0 . 

Recalling  the  result  cited  at  the  conclusion  of  §4,  we  may  here 
assume  that  all  the  A's  are  distinct,  and,  without  loss  of  gener- 
ality, take  X1  > X2  > A3  . Then  (5.2)2  yields  just  one  indepen- 
dent condition  which  can  be  expressed  in  the  form 

f"  = -Mf ' -e)/(X1-X3)2  , with  e > 0 . (5.4) 

The  substitution  from  (5.4)  into  (5.3)1  yields 

a = 4(A2n2  + ##) f ’ {Ff  ' + e[(X2-X2)2n2n2+  ## ] }/ (A^A^ 2 , (5.5) 

where 

F = (X1-X3)2(x2n2+  ##)  - {(A2-A2)2n2n2+  . (5.6) 

According  to  (5.5)  and  (5.6)  we  have  a > 0 if  and  only  if 

e > 0 provided  F >.  0 for  all  n and  F = 0 for  some  n . That 

these  latter  conditions  are  indeed  met  can  be  seen  by  expressing  F 
as  a function  of  n2  and  n2  (recall  n2=  l-n2-n2)  and  then 

looking  for  a minumum  of  F over  the  interior  of  the  triangular 

2 2 2 2 

domain  in  the  n^^-plane  bounded  by  the  lines  n^ 0,  n3=0, 

2 2 2,2 

n1+n3=l.  The  relevant  conditions  are  SF/Bn^  = 0 and  3F/3n2  = 0, 

which  lead  to  the  equations  [3] 
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(X1  "X2)ni  + CX3"X2)n3 


v 1 2J  1 ^3  2 3 


{Xl'X2'(Xl'X3)2}/2  ’ 


{X2'X3'(Xl'X3)2}/2 


(5.7) 


Since  X^  > X^  , these  are  found  to  be  inconsistent,  and,  hence, 
no  minimum  of  F exists.  It  is  also  found  that  F > 0 for 
n^  = 0 and  n^  = 0 , but,  for  n^+n^=l  , we  have 


F = (^1-x3)2(x1n1-  x3n3)3  » 


(5.8) 


which  vanishes  when 


ni  = X3/(X1+X3)  » n3  = X1/(W 


(5.9) 


We  conclude  that  F > 0 for  all  n except  for  the  four  values 
arising  from  (5.9),  for  which  F = 0 . 
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6 . Bifurcation  of  a Thick  Plate 

We  consider  a rectangular  plate  having  sides  of  length 
2 Z±,  2 Z2,  21^  and  which  occupies  the  region 

< ZA  < , A = 1,2,3  (6.1) 

in  its  undeformed  state.  The  plate  is  held  in  equilibrium  in  a 
state  characterized  by  (2.2)  by  means  of  normal  forces  uniformly 
distributed  over  the  faces  perpendicular  to  the  x1  and  x^ 
directions.  No  forces  are  applied  to  the  faces  initially  normal 
to  the  Xg-direction.  Here  we  investigate  the  possibility  that  a 
neighboring  equilibrium  state  can  exist  in  the  form  of  a small 
superposed  displacement  in  the  x1x2-plane  under  the  conditions  that 
no  further  deformation  occurs  in  the  x^ -direction , that  no  tangen- 
tial tractions  are  applied  to  the  load-carrying  faces,  and  that 
these  faces  remain  parallel  to  their  original  directions.  If  such 
bifurcation  is  possible,  we  say  that  the  underlying  pure  homogene- 
ous deformation  is  a critical  state. 

With  X ^ fixed,  we  note  from  (2.6)  and  (2.8)  that 

Ix  = A2  + (X  + A-1)**1,  I2  = A^X  + X"1)  + X"2  (6.2) 

where 

A = A2/X1  . (6.3) 


Also,  from  (2.14),  (2.20)  and  the  fact  that  the  faces  ?2  = ± 
are  force  free,  we  obtain 


and 
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irix  - 2(1-A 2)(W1  + X^2)/CXX3^  . (6.5) 

From  (4.1)^,  with  A * 3,  we  see  that  the  force  applied  to  the 
face  £ = 4 is  a thrust  or  tension  accordingly  as  A > 1 or 
A < 1 , respectively. 

To  describe  neighboring  equilibrium  states  we  take  u and  p 
in  the  form 

U1  = u1^1»e2^»  u2  = U2(Cl»^2)»  u3  = 0,  p = p(£1,£2)  . (6.6) 

Then  from  (2.18)  and  (2.9)  we  obtain 


J 


22  - 


We  employ  (2.16),  (6.3),  (6.4)  and  find  that  solutions  of 
(6.7),  (6.8),  (6.9)  can  be  obtained  in  the  separable  form  [4], 


-sin*??. 


cos^C- 


cos4>£. 


u!  " >Ui(£2),  u = }>U(£  ),  P = ^Q(52)  , 

cosH-J  1 2 sin’J'51  [ sin'i'C1  ‘ 2 


(6.10) 


where 


$ = niT/Jl1,  Y = (n-^Jir/J^,  n = 1,2,3,...  (6.11) 


and 


U1  " to  u’  * Q " ^2  {c12u"-n2(Cn+  e11-xe12)u-}  , 


(6.12) 


provided  that  U satisfies  the  differential  equation 


U^iv) -U2  + 1+A(X-1)2  }f22U"  + X2flUU  = 0 


(6.13) 


and  the  boundary  conditions 


U"  ♦ X2ft2U  = 0 


U"’  + S12(2X2  + 1+A(X-1)2}u'  = 0_ 


^2  " * *2 


(6.14) 


The  upper  set  of  expressions  in  (6.10)  represents  a symmetric  mode 
of  deformation  with  respect  to  the  x^x^-plane  and  the  lower  set 
an  anti-symmetric  mode  with  respect  to  this  plane.  Also,  the  quan- 
tity n denotes  either  $ or  'V  , depending  on  which  set  is 
considered,  and  A is  defined  by 


(6.15) 


It  appears  that  Wesolowski  [14]  was  the  first  to  exploit  the 


J 
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fact  that  the  governing  differential  equation  and  boundary  con- 
ditions for  a bifurcation  problem  of  this  type  depend  on  the  strain- 
energy  function  only  through  a single  material  parameter,  such  as 
A . We  remark  here  that,  for  any  specified  value  of  , A is  a 
function  of  X only  (cf.(6.2)}  . 

With  (3.15)  we  see  that  A can  be  expressed  as 

A = (X+l) 2M^/ (XX^K^)  and,  from  (4.1)  (with  the  subscript  A=3)  , 
together  with  (2.6),  (2.8)  and  (6.3),  that  the  Hadamard  criterion 
can  be  satisfied  only  if  (X-l)  2M^/  (XX^^)  > -1  . Whence,  A must 
satisfy  (X-1)2A  > -(X+l)2  , or  simply 

A > - (X  + l) 2/  (X-l) 2 = Aq(X)  , say  , (6.16) 

in  order  that  the  underlying  deformation  meet  the  minimal  stability 
requirement.  Of  course,  no  restriction  arises  from  (6.16)  if  X=1  . 

Upon  assuming  a solution  of  (6.13)  of  the  form  U = exp(ftr£2)  , 
we  obtain  a biquadratic  equation  for  T which  has  four  roots 
of  the  form  ± (R±T)  where 

R = 3s(  (X  + l) 2 + A(X-l)2}is  , T = (A+l)  (X-l)2}^  . (6.17) 


According  to  (6.16),  R is  real  for  all  possible  values  of  A and 


T is  real  if  A >.  - 1 , 

T is  imaginary  if  -1  > A > - (X  + l) 2/ (X-l) 2 . 


(6.18) 


It  is  convenient  to  treat  these  two  cases  separately.  The  first 
is  discussed  in  § 7 and  the  other  in  §8. 
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7 . Bifurcation  When  A > - 1 

When  (6.18)^  applies,  we  let 


r1=R+T,  t2=r-t 

and  write  the  general  solution  of  (6.13)  in  the  form  [4] 


U = L1coshnriC2  + LgSinhflr  £2 

+ N1cosh^T2C2  + N2sinhar2£2  , ( 

where  L^LgjN  ,N  are  constants.  The  substitution  into  (6.14) 
yields  four  equations  for  the  determination  of  these  constants, 
which,  with  the  notation 


n = a*2  , 


can  be  written  as 


(r^+x2)L1coshr1n  + (r2+x2)N1coshr2ri  = o , 
r1(r2+x2)L1sinhr1n  + r2(r2+x2)NlSinhr2n  = o , 
(r2+x2)L2sinhrln  + (r2+x2)N2sinhr2n  - o , 
rl(r2+x2)L2coshr1n  + r2(r2+x2)N2coshr2n  = o . 


In  order  that  non-trivial  solutions  of  (7.4)  exist,  we  must 


have  either 


L„  » 

2 2 


= 0 and 


tanhF^n 

tanhr, n 


r|n^  \ 2 
r?*x2  , 


Li  * Ni 


= 0 and 


tanhf^n 

tanhr2n 


r2+x2  \ 2 


_1  

r2  r2n2 
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If  (7.5)  holds,  we  can  obtain  non- trivial  values  of  L1  and  N1 
from  (7.4)  and,  accordingly,  from  (6.10),  (6.12)  and  (7.2)  we  see 
that  the  variation  of  U-j^  is  odd  and  u2  is  even  with  respect  to 
£2  . This  gives  rise  to  an  anti -symmetric  shape  with  respect  to 
the  x^x^-plane,  which  we  call  a flexural  mode.  Conversely,  the 
shape  associated  with  the  satisfaction  of  (7.6)  is  symmetric  with 
respect  to  the  x^^-plane  and  we  call  that  a barreling  mode. 

With  the  convention  that 


v 


for  a flexural  mode 
for  a barreling  mode  , 


(7.7) 


and  with  (7.1),  it  is  possible  to  write  both  bifurcation  conditions 
(7.5)2  and  (7. 6) 2 in  the  form  [4] 

sinh(2Rri)  _ R { 4T2  + X (A  + l) 2 } ,,  0, 

■ V T UR2.Hi.^}  ■ (7-8) 

From  (6.17),  and  bearing  in  mind  that  A is  a function  of  X , we 
see  that  R and  T are  functions  of  X . Then  (7.8)  can  be  re- 
garded as  an  equation  for  the  determination  of  n as  a function  of 
X . From  (6.11)  and  (7.3), 

n = n-ni^/21  , n = 1,2,3,...  (7.9) 

where  n denotes  the  number  of  half -wavelengths  parallel  to  the 
xi-direction  in  a neighboring  mode.  Thus,  we  may  also  regard  (7.8) 
as  an  equation  for  the  determination  of  critical  values  of  X 
corresponding  to  a particular  aspect  ratio  ^2^1^  anc*  a §^ven 
mode  shape. 
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ain  results  based  on  (7.8)  have  been  proven  analytically 
the  assumption  that  A is  constant . These  are  summar- 
ollows,  where  9 = 4R2-X(X-1)2  : 


If  X > 1 (tensile  case)  (7.8)  admits  no  real 
solutions  and,  hence,  no  bifurcation  is  possible. 

For  any  X > 1 , one  and  only  one  positive  value  of 
n exists  for  flexure  if  9 > 0 and  none  exists  if 
9 < 0 . 

For  any  X > 1 , one  and  only  one  positive  value  of 
p exists  for  barreling  if  9 < 0 and  none  exists  if 
9 > 0 . 


2 2 

the  vanishing  of  (4R  -X(X-l)  } is  the  condition  which 
the  ranges  of  critical  values  of  X for  which  flexural 
ling  bifurcations  can  occur.  From  (6.17)^,  this  separa- 
tion is  defined  by 


A = X-(X+1)2/(X-1)2  = A (X) , say  . 


(7.10) 


be  said  about  this  in  the  next  section.  The  function  A 
plotted  in  the  XA-plane  in  Fig. 3 of  [4]  for  A > -1 


rical  results  from  (7.8)  can  be  obtained  [4]  by  assigning 
on.3tant  values  to  A and  then  solving  numerically  for  n 
tion  of  X . This  procedure  yields  curves  of  the  form 
ematically  in  Fig.l.  here  we  have  taken  three  constant 
,b,c  for  A with  the  assumed  ordering  -1  < a < b < c . 
curves  represent  solutions  of  (7.8)  for  flexure  and  the 
rves  the  solutions  for  barreling. 

ough  the  condition  A = constant  is  somewhat  artificial 
of  course,  for  the  Mooney-Rivlin  material  for  which  A = 0) , 
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from  curves  such,  as  those  in  Fig.l  it  is  possible  to  determine 
critical  values  of  X , for  any  specified  value  of  n , for 
arbitrary  variation  of  A with  X . The  procedure  is  to  read 
off  pairs  of  values  of  critical  X and  corresponding  A for  the 
value  of  n specified  and  then  plot  these  points  in  the  XA-plane. 
This  yields  two  distinct  curves,  one  for  flexure  and  one  for  bar- 
reling. On  this  same  graph  we  plot  the  actual  variation  of  A 
with  X from  (6.15)  (recall  (6.1)  and  the  fact  that  X^  is  fixed}. 
The  intersections  of  this  latter  curve  with  the  former  give  the 
critical  values  of  X for  flexure  and  barreling  corresponding  to 
the  value  of  n for  which  the  original  curves  were  derived.  This 
procedure  is  illustrated  in  Fig. 4 of  [4]. 

We  take  this  opportunity  to  correct  equation  (7.4)  of  [4] 
which,  for  the  example  chosen,  should  read 

A = 2(X+1)2/(X2+6X+1)  . (7.11) 

This  correction  alters  slightly  the  numerical  values  stated  in  the 
final  sentence  of  [4] , but  does  not  change  in  any  way  the  qualita- 
tive aspects  of  the  discussion  there.  Of  course,  if  the  actual 
variation  of  A with  X is  known,  as  for  example  in  (7.11),  it 
could  be  employed  in  (7.8)  and  a numerical  solution  effected 
directly . 
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8 . Bifurcation  when  A < - 1 

When  (6.18)~  applies,  we  write 

T = iS  , S = 3s(-(A+l)  U-l)2}1"2  , (8.1) 

and  express  the  general  solution  of  (6.13)  in  the  form 


U = (I^cosflSCg  + L2sinf2SC2)  cosh^R52 

+ (N1sinJQS?2  + N2cosnS?2)  sinh^R?2  , 


(8.2) 


where  Li>L2,N1,N2  are  constants  and  R is  given  by  (6.17)^. 
We  introduce  the  notation 


q-L  ■ R2-  S2  + A2;  q2  = 2RS ; q3  = (l-X)R;  q^  = (l-A)S  ; 

Y = coshRn  cosSn  ; Z = sinhRn  sinSn  ; (8.3) 

y = coshRn  sinSn  ; z = sinhRn  cosSn  ; 


where  n is  given  by  (7.3),  and  then  substitute  from  (8.3)  into 
(6.14).  This  yields  {cf.(7.4)> 


L1(qiY-q2Z)  + N1(q1Z+q2Y) 

Mq3Z  + qUy)  + N2(q3y-quz) 

^2(<\1y  + ^2Z'>  + N2CqiZ_q2y) 

L2(q3Z-quY)  + N2(q3Y+quZ) 


(8.4) 


(8.5) 


Now  (8.4)  and  (8.5)  admit  non-trivial  solutions  for  the  L's 
and  N's  only  if  either 

L2  = N2  = 0 and  (q1Y-q2Z)  (q3y-q1|t)  = (q^+qgY)  (q3z  + quy)  (8.6) 

or 

... ;;  ~~  ■ • — 
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L]_  = Nx  = 0 and  (q^+q^)  Cq3Y+quZ)  = (q^ -q2y)  (q3Z -q^Y)  , (8.7) 

unless  it  should  happen  that  the  determinants  associated  with 
(8.4)  and  (8.5)  both  vanish  simultaneously;  that  is,  unless 


(Yz+Zy) (q1qu+q2q3)  = 0 . 


(8.8) 


The  significance  of  this  exceptional  case  shall  become  apparent 
below. 

If  (8.6)  applies  then  non-trivial  values  of  L1  and  N1  can 
be tfound  from  (8.4),  which,  with  (6.10),  (6.12)1  and  (8.2),  gives 
values  for  u1  and  u2  whose  variations  with  £2  are,  respec- 
tively, odd  and  even.  As  before,  this  situation  describes  flexural 
modes.  The  converse  behavior  is  noted  when  (8.7)  applies  which, 
therefore,  governs  barreling  modes. 

We  let 


r = X (3-X)  (1  + X)  + (A+l) (X-l)2  , 
s = X (X  + l)  + (A+l)  (X-l)2 


(8.9) 


and  substitute  from  (6.17),  (8.1)  and  (8.3)  into  (8.6)  and  (8.7) 
and  cast  both  bifurcation  conditions  in  the  form 


vRs  sin  2Sn  = Sr  sinh  2Rn  , (8.10) 

where  v is  given  by  (7.7).  Bearing  in  mind  that  R,S,r,s  are 
functions  of  X , we  have,  as  with  (7.8),  that  (8.10)  serves  to 
define  a relationship  between  n and  critical  values  of  X . 

It  can  readily  be  shown  that  (8.8)  is  satisfied  trivially  if 
X » 1 or  if  n = 0 ; but,  for  X ^ 1 , there  are  non-trivial 
solutions  for  n of  the  form 
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n = kn  Ck=l,2 ,3, . . . ) , fj  = tt/2S  , (8.11) 

if  and  only  if 

r = 0 . (8.12) 

The  solutions  (8.11)  also  satisfy  (8.10)  and,  with  (8.12),  are 
seen  to  be  independent  of  the  value  assigned  to  v . This  leads 
to  the  possibility  of  the  simultaneous  appearance  of  flexural  and 
barreling  modes,  a feature  which  stands  in  sharp  contrast  to  that 
when  A > - 1 , where  these  modes  are  clearly  delineated.  The 
condition  (8.12),  which  gives  rise  to  this  feature,  is  seen,  with 
(8.9)^,  to  yield  an  expression  for  A in  terms  of  X that  is 
identical  in  form  to  A in  (7.10). 

Bearing  in  mind  that  X^  is  fixed,  we  employ  (6.2),  (6.5), 

(6.15)  and  calculate 

dnxi/dX  = -(W1  + x2w2){3X2  + l+A(X-l)2}/(X3X3)i£  , (8.13) 

and  note  that  II 7 n is  stationary  for  any  value  of  X which 
satisfies 

A = - (3X2+1) / (X -1) 2 = A(X)  , say  . (8.14) 


From  (6.16),  (7.10),  (8.14)  it  easily  follows  that 


> 

o 

A 

> 

< A 

for 

T~i 

V 

V 

o 

A < Aq 

< A 

for 

X > 1 

(8.15) 


Plots  of  AQ , A,  A in  the  XA-plane  are  shown  in  Fig. 2.  The 
region  below  the  curve  A^  represents  values  of  A which  violate 

(6.16)  and,  thus,  is  irrelevant. 
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To  proceed  with  the  discussion  of  (8.10)  we  again  resort  to 
the  somewhat  artificial  situation  where  A is  assigned  some  con- 
stant value  < - 1 . In  doing  so,  from  Fig. 2,  we  see  that  the 
possible  values  of  A must  be  restricted  to  lie  in  the  interval 
e £ X £ E where  e and  E are  the  roots  of  the  equation 
A = Aq(x)  . That  is,  with  (6.16), 

e = -(1-/^A)2/(A+1)  , E = -(1+/^A)2/(A+1)  . (8.16) 

We  first  consider  the  case  where  n is  a thrust  (A>1)  and  then 
where  it  is  tensile  (A<1)  . 

8.1  The  Case  of  Thrust 

Here,  for  a specified  value  of  A < -1  , we  consider  values 
of  X that  lie  in  the  interval  1 £ X £ E . Let  X denote  the 
abscissa  of  the  intersection  of  the  horizontal  line  corresponding 
to  this  value  of  A with  the  curve  A in  Fig.  2.  Then  1 < X < E , 
and  we  note  that  X and  E both  increase  monotonically  with  in- 
creasing A . 

The  solutions  of  (8.10)  that  result  by  assigning  two  values 
to  A are  shown  schematically  in  Fig.  3,  where  values  of  n are 
plotted  against  critical  values  of  X . The  solid  curves  pertain 
to  flexural  bifurcation  and  the  broken  curves  to  barreling. 

The  numerical  values  of  the  relevant  parameters  X,  E and  n 
corresponding  to  three  representative  values  of  A are  listed  in 
Table  I.  These  can  be  used  to  gain  some  idea  of  the  scale 
associated  with  Fig. 3.  It  can  be  shown  that  the  minimum  value  of 
n is  approximately  2 tt/ 5 . 
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Table  I.  Relevant  parameters  for  three  values  of 
A (thrust) 


A 

X 

E 

n 

— 

-2.4625 

2.6 

4.5135 

1.6236 

-7.0 

2.0 

2.2153 

1.2825 

-17.7778 

1.6 

1.6360 

1.3018 

The  features  common  to  the  solutions  of  (8.10),  for  any  con- 
stant A < -1  , are  summarized  as  follows.  Each  flexural  curve 
emanates  from  the  point  X = 1 , n = 0 and  increases  steadily  to 
a point  X = X , n = n where  it  intersects  a downward  sloping 
barreling  curve.  This  barreling  curve  decreases  steadily  until  it 
intersects  a limiting  line  X = E at  a value  of  n,  nQ  , say, 
where  T > nQ  > 0 . No  barreling  bifurcation  is  possible  if 
n < nQ  • The  flexure  and  barreling  curves  intersect  at  the  points 
(cf.(8.11)}  X = X , n = kn(k=l,2 ,3, . . .)  , and  each  curve  crosses 

the  line  X = X at  these  points.  If  n lies  in  any  of  the  open 

intervals  (IT,  2rf)  , (3n,4'n),  (5rT,6n),  etc.,  the  critical  values  for 
barreling  are  less  than  X , while  those  for  flexure,  to  the  extent 
they  exist,  are  greater  than  X . Conversely,  for  values  of  n in 
the  intervals  (0,n),  (2n,3n),  (4n,  5T),  etc., the  flexure  curve  lies 
to  the  left  of  the  line  X = X , and  the  barreling  curve  to  the 
right.  Both  the  flexure  and  barreling  curves  have  the  line  X = X 
as  an  asymptote  as  n ■*  ® . The  barreling  and  flexure  curves  may 

intersect  the  line  X = E ; the  tendency  for  such  intersections  to 

occur  increases  as  A decreases  (algebraically) . 

Recalling  that  n is  given  by  (7.9),  it  appears  that  certain 
peculiarities  arise  in  connection  with  the  interpretation  of  the 
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solutions  depicted  in  Fig. 3 for  values  of  X "close"  to  X . If 
taken  at  face  value,  these  results  for  large  n can  be  considered 
as  conditions  for  the  appearance  of  surface  modes  in  a thick  plate 
[10].  However,  pending  further  investigation  into  the  relevance  of 
the  critical  conditions  for  such  plates,  we  here  eschew  giving  any 
direct  physical  interpretation. 

On  the  other  hand,  for  values  of  X near  unity,  the  flexural 
curves  in  Figs.l  and  3 show  a strong  similarity.  This  region, 
which  gives  the  bifurcation  condition  for  a thin  plate,  has  been 
studied  [5]  for  the  case  A > -1  , governed  by  (7.8).  The  result 
is  expressed  by  the  asymptotic  series 

X = 1 + a;Ln2+  a2T)k+  a3n6+  aun8  + 0(n10)  , (8.17) 

where  the  a's  are  given  by 


ax  = 2/3  , a2  = 16/45  , a3  = 2 (A  +18/7)/27  , 
a^  = 2(A1+113/175)/27  , 


(8.18) 


with  A1  denoting  the  value  of  A when  X = 1 . From  (8.10), 
following  the  expansion  procedure  used  in  [5],  we  again  obtain 
(8.17),  with  (8.18).  This  result  shows  (i)  that  the  critical  value 
of  X for  flexure  is  independent  of  both  the  form  of  the  strain- 
energy  function  and  the  initial  extension  ratio  X^  up  to  the 
terms  of  fourth  order  in  n and  (ii)  that  this  critical  value  is 
extremely  insensitive  to  the  detailed  behavior  of  the  material  up 
to  terms  of  eighth  order  in  n . 

Numerical  results  for  arbitrary  variation  of  A with  X 
(i.e.,  for  an  arbitrary  material}  could  be  obtained  by  plotting 
curves  as  in  Fig. 3 for  a large  number  of  constant  values  of  A 
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and  then  following  the  procedure  described  in  [4]  . Alternatively, 
the  function  A could  be  employed  directly  in  (8.10). 

8.2  The  Case  of  Tensile  Load 

Here,  for  a specified  value  of  A < -1  , we  consider  values 
of  X that  lie  in  the  interval  e < X < 1 (cf.(8.16)}  , and  let 
X,X  denote  the  abscissae  of  the  intersections  of  the  horizontal 
line  corresponding  to  this  value  of  A with  the  curves  A,  A in 
Fig. 2.  Then  e < X < X < 1 , and  we  note  that  e,  X,  X all 
decrease  monotonically  as  A increases. 

In  Fig. 4 are  shown,  schematically,  the  solutions  of  (8.10) 
that  result  by  assigning  two  values  to  A . Again,  the  solid 
curves  pertain  to  flexure  and  the  broken  curves  to  barreling. 
Wesolowski  [14]  restricted  attention  to  barreling  bifurcations. 

To  obtain  some  idea  of  the  scale  associated  with  Fig. 4,  the  numeri- 
cal values  of  e,  X,X,n  , corresponding  to  three  representative 
values  of  A are  listed  in  Table  II. 

The  most  important  feature  to  note  in  connection  with  any 
solution  of  the  type  shown  in  Fig. 4 is  that  neither  the  barreling 
nor  the  flexural  curve  has  any  point  lying  to  the  right  of  the  line 
X = X . Thus,  there  is  a range  of  X , X < X < 1 , for  which  no 
bifurcation  can  occur. 


Table  II.  Relevant  parameters  for  three  values  of 
A (tensile) 




A 

e 

A 

X 

n 

-2.05 

0.1776 

0.2 

0.2425 

3.8324 

-7.0 

0.4514 

0. 4641 

0 . 5 

2.3933 

-15.4 

0.5938 

0.6 

0.6246 

2.0697 

If  we  assume  for  the  moment  a material  does  exist  for  which 


k 


J 
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A is  constant  and  which,  also  meets  the  minimal  material  stability 
requirement  (W^X^)  > 0 > then  from  (6.5),  (8.13)  and  (8.14) 
it  follows  that  the  tensile  force  II  increases  from  zero  to  a 
maximum  as  X decreases  from  unity  to  X , and  then  decreases  as 
X continues  to  decrease.  Whence,  in  concurrence  with  Wesolowski 
[14]  and  Hill  and  Hutchinson  [9] , we  conclude  that  if  a bifurcation 
is  to  occur,  it  must  do  so  after  a load  maximum  is  achieved.  How- 
ever, this  behavior  can  be  used  to  illustrate  that  a plate  of  this 
material  having  a finite  aspect  ratio  (J^/J^)  could  behave  as  a 
perpetual  motion  machine  without  suffering  a bifurcation. 

The  extension  of  this  argument  to  a general  material  leads  to 
the  conclusion  that,  if  the  relevant  variation  of  A with  X is 
plotted  in  Fig.  2,  the  resulting  curve  must  lie  above  the  curve  A 
for  0 < X < 1 . Thus,  in  general,  no  bifurcations  of  the  types 
considered  can  exist  under  tensile  loading  conditions. 

We  remark  that  the  curve  A , itself,  is  independent  of  the 
value  of  the  extension  ratio  X^  . Also,  the  above  conclusion  is 
based  solely  on  the  implications  arising  from  (6.5),  (8.13)  and  the 
necessity  of  the  Baker-Ericksen  inequalities,  and  leads  to  a further 
restriction  on  the  form  of  the  strain-energy  function  beyond  that 
implied  by  (4.1)2*  This  restriction,  viz.  A > A , can  be  cast 
in  the  form 

(X2-X2)2 M3  + (X2+3X2)K3  > 0 , (X1>X2)  (8.19) 


upon  employing  (2.8),  (3.15),  (6.3),  (6.13)  and  (8.14).  In  the 
event  that  the  passive  load  applied  to  the  faces  £3  = ± A 3 
of  the  plate  is  also  tensile  (i.e.,  if  X3>X2)  , then  a similar 
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the  roles  of  X.^  and  X^  interchanged, 
that  another  of  the  conditions  (4.1)? 
stronger  restriction 

- (*2+3x2)!^  > 0 , (X3>X2)  . (5.20) 
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Plots  of  n vs.  critical  X for  two  typical  values  of 
under  tensile  loading. 
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